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Abstract 

We investigate the feasibility of explosive particle production via parametric resonance or tachy- 
onic preheating in multi- field inflationary models by means of lattice simulations. We observe a 
strong suppression of resonances in the presence of four-leg interactions between the infiaton fields 
and a scalar matter field, leading to insufficient preheating when more than two inflatons couple 
to the same matter field. This suppression is caused by a dephasing of the inflatons that increases 
the effective mass of the matter field. 

Including three-leg interactions leads to tachyonic preheating, which is not suppressed by an 
increase in the number of fields. If four-leg interactions are sub-dominant, we observe a slight 
enhancement of tachyonic preheating. Thus, in order for preheating after multi- field inflation to be 
efficient, one needs to ensure that three- leg interactions are present. If no tachyonic contributions 
exist, we expect the old theory of reheating to be applicable. 
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I. INTRODUCTION 



After inflation, the universe needs to heat up to temperatures above lMeV in order 
for primordial nucleosynthesis to commence. This temperature increase could occur via 
the perturbative decay of the inflaton(s) as described by the old theory of reheating [l[ 
2], or by the rapid production of particles during a phase of preheating (see |3j, |4| for 
reviews). The latter mechanism usually involves resonances or instabilities, with a resulting 
temperature well above lMeV, sometimes as high as 10 15 GeV, as opposed to the moderate 
temperatures accessible by means of the old theory of reheating 94J. High temperatures may 
be advantageous for scenarios of GUT baryogenesis (see i.e. |5J), but are usually in conflict 
with relic-bounds, one of which originates, for instance, from a prevention of gravitino- 
overpro duct ion (T r h < 10 9 GeV). For these reasons it is crucial to understand whether 
or not preheating occurs after inflation. 

Our current knowledge of preheating is primarily based on field theoretical models describ- 
ing the inflatons' decay by non-perturbative effects, such as parametric resonance [s, [9], [h]] 
or tachyonic preheating [ll] among others (see jjj, 4] for reviews). This approach should be 
seen as a preliminary step in the absence of a better understanding of the degrees of freedom 
that drive and end inflation. For instance, it may be that stringy effects are crucial, such 



as in brane inflation models where inflation ends once branes annihilate [14[. The hope is 
that an effective field theoretical approach emerges, once we have a better understanding of 
the microphysical nature of inflation. In the meantime, it is desirable to have a thorough 
understanding of the conditions under which the energy transfer from the inflationary sector 
to preheat matter fields is efficient. Further, the preliminary theory of preheating has been 
successfully applied in other areas, for instance in the notion of moduli trapping on the 



landscape 



15 



16 



171 ]. a possible solution [18| to the moduli-induced gravitino problem [95 ] 



or as an additional source of gravitational waves 
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26j. In light of 



these applications, it is necessary to investigate the viability of the proposed mechanisms, 
particularly as, in recent years, new developments in string cosmology (see 27, 28, Q, [30} 
for recent reviews) have led to the emergence of several multi-field inflationary models; a 



33, 34] 
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see also 
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ew examples are inflation on the landscape 3l|, |32j, inflation driven by axions i.e. A/"-flation 
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36] ), inflation driven by tachyons 37], |38| or by multiple M5-branes 



421 ] among others (see 43J for a review). The majority of these models are 
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realizations of assisted inflation 



44 



45 



46l . 1471 . |48| |. where many fields are dynamical. An 



advantage of these models is the possible alleviation of the 77-problem, since any given field 
does not need to traverse a large distance in field space, given that many fields (J\f ~ 10 3 ) 
assist each other in driving the inflationary phase. 

Unfortunately, little is known about the actual process of reheating standard model (SM) 
degrees of freedom in almost all inflationary scenarios, and there is a danger of primarily 
heating hidden sectors, especially if many fields contribute to inflation (see 49j for a case 
study within AA-flation). If inflation would be driven with the field content of the MSSM 



50. 
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preheating 



52 



the standard model can be properly reheated 52j, for instance via instant 



541 ]. since the couplings of these fields are known. However, it is challenging 



to achieve the needed sixty e-folds of inflation within this scenario; as a result, a prior phase 
of high-scale inflation is usually assumed, commonly terminated by a phase of preheating. A 
related possibility to reheat primarily standard model particles, which does not involve low 
scale inflation, consists of a second phase of reheatin g, ca used by the decay of a long lived, 
light particle species such as the saxion (see i.e. [55|, |56| 96]), which dilutes any unwanted 
relics and overabundance of hidden sector particles produced during preheating. 

However, even a field theoretical description of preheating is by no means complete; we 
lack an understanding to which fields the inflatons couple, what type of couplings are present, 
and how big the coupling constants are. Non-perturbative mechanisms are common, given a 
certain amount of fine tuning, but they have mainly been studied in simple setups involving 



one 



EH, 



63 



64J or two 



65. 



66 1 inflatons coupled to a single matter field, or one inflaton 



671] (see however the recent case study of preheating after 



coupled to several matter fields 
A/"-flation with J\f ~ 100 fields [68|, |69j). Preheating is complicated by the importance of 
backreaction after the first few bursts of particle production which affects resonances; these 
effects can only be recovered using lattice simulations 



70 



71 



,172J,|73|,|74J, which have only 



become feasible for more complicated scenarios in recent years. Further, even if the inflatons 
decay via this route, the universe is still in a non-thermal state. The era of thermalization 



is involved and again requires lattice simulations [75l. 1 76]. 



In this article, we investigate the early stages of preheating in multi-field inflationary 
models via a lattice simulation, incorporating backreaction and up to five inflatons. We 
use an implementation of LatticeEasy [70] (Defrost, an alternative code developed by 
A. Frolov, can be found in 72[]). Our aim is to investigate conditions under which para- 
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metric resonance (quadratic couplings between matter and inflaton fields, that is four-leg 
interactions) and tachyonic preheating (Yukawa couplings, that is three-leg interactions) 
are efficient. By efficient we mean that the majority of the energy in the inflaton fields is 
transferred to preheat matter fields within a few hundred oscillations. Our primary concern 
is the alteration of resonance effects or instabilities if more than one inflaton field is coupled 
to the same matter field. For instance, in the presence of a few fields, the expectation is 
that resonances are enhanced, since stability bands are generically destroyed, leading to 
instabilities of almost all Fourier modes of the matter field; this is sometimes referred to 
as Cantor Prehea ting 66j and has been studied numerically in a two-field model (without 
backreaction) in 65j. However, in the presence of many fields, inflatons generically run 
out of phase, which in turn increases the effective mass of the matter field so that efficient 
preheating becomes impossible |68|, |69(. Our first goal is to understand quantitatively the 



competition of these two opposing effects as we increase the number of fields one by one. 
Secondly, we examine tachyonic preheating (three- leg interactions), which is not expected 
to be hindered by dephasing effects 7?J ■ In the presence of both, three- and four-leg inter- 
actions, it is unclear if an increase of the matter fields effective mass leads to a suppression 
of instabilities. 

We find that resonances are strongly suppressed in the presence of four-leg interactions 
if more than two inflatons couple to the same preheat matter field; further, even in the 
presence of only two inflatons, resonances are not enhanced, as indicated by Bassett et.al. 



in 



65], but slightly suppressed. Thus, Cantor preheating is less efficient than previously 



anticipated. Turning our attention to Yukawa couplings, we find that particle production 
remains efficient. There is no suppression if the number of inflatons is increased, and if 
four-leg interactions are suppressed, tachyonic preheating is sightly enhanced. 

The outline of this paper is as follows: after introducing a generic but simple multi- 
field inflationary model in Sec. [Ill we provide the results of lattice simulations for standard 
parametric resonance models based on four-leg interactions in Sec. Ill Al mixed three- and 
four-leg interactions in Sec JIIBI and tachyonic preheating (three-leg interactions dominate) 
in Sec. Ill CI We conclude in Sec JIIII Details about our implementation of LatticeEasy 
can be found in Appendix [A] Throughout this article we set the Planck mass equal to one, 
mJ = l/G = l. 
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II. PREHEATING 



Consider several inflaton fields tfi, % = 1...M, with canonical kinetic terms. After 
inflation, each inflaton approaches its respective minimum in the potential. Since we are 
interested in preheating, we expand the potential around the minima so that, after redefining 
the fields such that the global minimum lies at = 0, the total potential becomes W = ^ Vi 
with Vi = mfipj/2] generically m 8 7^ rrij for i 7^ j, so that the fields oscillate with different 
frequencies. In simple models of inflation the slope during prehea ting is related to the slope 



during inflation. These models are also the most predictive ones 



78 



79 



801 ] . since more ad 



hoc features in the potential allow the fitting of more or less any exotic data. For instance, if 



inflation is driven by a single field with a quadratic potential (the simplest mode 
its mass is determined by the COBE normalization = (2.41 ± 0.11) x 10~ 9 



imaginable), 



8l| ; for more 



than one field with simple quadratic potentials the effective single field still needs to sa tisfy 



the COBE normalization, but the actual spread of masses is not significantly constrained [97]. 
Heavy fields become irrelevant during inflation, since they quickly roll down the potential. 
As a consequence, we expect a narrow spread of masses for fields relevant for preheating, at 
least in simple models (this is indeed the case in A/'-flation 68]). Based on this reasoning, 
we distribute mf equidistantly over the interval 0.5 x 10~ 12 ... 1.5 x 10~ 12 , i.e. m\ = 10~ 12 
for Af = 1, ml = 0.5 x 10" 12 and m\ = 1.5 x 10~ 12 for Af = 2, m\ = 0.5 x 10~ 12 , m 2 = 10~ 12 
and m\ = 1.5 x 10~ 12 for Af = 3 etc. Note that we ordered the masses: m ; < rrij ifi< j. We 
further define the average square mass as m 2 = 10~ 12 and a dimensionless time via r = tm 
(see Appendix [S} ; the latter one is used in all plots. 



In realistic scenarios, it may be that several fields decay during inflation 
assist inflation only during early stages [68 



31 



82 



83| or 



Such fields do not make a considerable con- 
tribution to preheating and are therefore excluded. The fields we consider should have 
comparable energies when preheating commences. For simplicity, we impose equal en- 
ergy initial conditions m 2 <£> 2 (0) = m 2 <£> 2 (0) and set the total potential energy equal to 
^m 2 ^ 2 (0)/2 = (10^ 6 ) 2 x (0.193) 2 /2 w 1.8 x 10~ 14 at the onset of preheating, the usual 
value for a single inflaton field with a quadratic potential [10]. These initial conditions 
simplify comparisons of setups with a varying number of fields. 

To model preheating, we couple all inflaton fields to the same scalar matter field x- m 
a realistic scenario we expect a more complicated matter sector, but focusing on only one 
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matter field should provide an instructive toy model. We allow for Yukawa type interac- 
tions aiifiX 2 /2 and quadratic interactions gitp 2 x 2 /2, since both are generically present [98| ; 
however, we ignore terms higher order in (pi as well as cross coupling terms proportional 
to (fi^Pj for i j (the former ones are sub-leading in the Taylor expansion and the latter 
ones are usually suppressed for multi-field inflationary models of interest, such as assisted 



inflation 



44|). During preheating, only light degrees of freedom are produced [7l|, |84|, hence, 



we neglect the necessarily small bare mass of \. However, to guarantee that the potential 
is bounded from below for large x an d to increase the stability of our code, we include a 
small self interaction term A% 4 /4; this term prevents an unphysical runaway behavior of \ 



in our numerical integration without influencing preheating significant 



coupling constants between the inflatons and x to be comparable 
Oi = o. All in all, the potential we consider during preheating reads 



y. Because we expect 



99| | . we take ^ = g and 



x 



(i) 



The equations of motion are the Friedmann equations which we combine to 



d + 2- 



«7T 

a 



(2) 



as well as the Klein-Gordon equations 



(pi + 3-tpi 
a 



a 1 n 

X + 3-x-^V 2 x + 
a nZ 



dW 

d<pi 
dW 



0. 

= 0. 



a" dx 

We integrate the above set numerically using the LatticeEasy code 



(3) 
(4) 



701 ] which uses a 

staggered leapfrog integrator. The computation is performed in real space on a n 3 = 128 3 
lattice, with a comoving box size of L = 5/m, while the initial conditions (including small 
initial inhomogeneities) are prescribed in Appendix |X] (see also This means we can 

cover wave-numbers in the interval 2w/L < k < 2-n v3n / '(2L) . 

In the next sections, we discuss the results of various choices of M and the coupling 
constants a and g. For more details on our implementation of LatticeEasy see Appendix 

m 



7 



Case 


g x (0.193) 2 /m 2 


cr x 0.193/m 2 


A x (0.193) 2 /m 2 


A 


10 4 





5 x 10 3 


B 


10 4 


100 


5 x 10 3 


C 


100 


100 


10 4 



TABLE I: The coupling constants for cases A (four- leg interactions, Sec. Ill AT) . B (three- and four- 
leg interactions, Sec. Ill Bp and C (three-leg interactions and minor four-leg interactions (added for 
stability), Sec. Ill Cj) . A ^ suppresses an unphysical runaway behavior of \i but does not interfere 
with preheating; here m = 10 -6 and the Planck mass is set to one. For ease of comparison we 
chose values identical to those in [ll|], where tachyonic preheating for M = 1 is discussed (note 



that instead of g 2 as in 



111 ] we use g in the potential). 



A. Parametric Resonance (a = 0, g / 0): Suppression due to Dephasing 



In the absence of three-leg interactions and considering a single inflaton fie 



d (M 



we 

observe the well known amplification of x due to parametric resonance 0, H hi (Fig. m 
A and [2] A, M = 1). To be concrete, if backreaction and expansion effects are ignored, 
the equation of motion for the Fourier modes of the matter field Xk turns into the Mathieu 
equation 



x'l+(A k -2q C os(2r)) X k = 0. 



(5) 



where a prime denotes a derivative with respect to r = mt, Ak = k 2 jm 2 + 2q, q = g$ 2 / (4m 2 ) 
and $ is the initial amplitude of the inflaton field oscillations (p = $cos(r); the Mathieu 
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851 ]. resulting in narrow 



equation exhibits well known stability and instability bands 

or broad resonance [h]]. If a mode lies within an instability band, its amplitude grows 
exponentially cx e^ k , where Hk > is the Floquet index. Including the expansion of the 



universe causes modes to shift through bands, leading to stochastic resonance [lOj. Here, 
one can still derive approximate analytic expressions describing the growth of the matter 
field [lOj, given that only one inflaton field is involved. However, in all cases backreaction 
becomes important soon after the first spurts of particle production, causing the inflaton 
field to fragment and leading to an extended turbulent regime followed by an even 

longer phase of thermalization Q, Q| • We checked that we recover these known results in 



our simulations (see for instance Fig. [2] for the variance (TA22|) in the N = 1 case). 
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In the presence of two fie 



on spectral theory 



86. 



S3, 



88 



ds, Bassett argued that resonances could be enhanced: based 



89 



90 



91 



921 ]. he showed that the stability bands generically 



dissolve into a nowhere dense set, similar to a Cantor set, if the ratio of the inflatons' 
oscillation frequencies is not a simple fraction; then almost all modes should be amplified 



66f| . This is sometimes referred to as Cantor preheating. However, no analytic expressions 



are known for the magnitude of the generalized Floquet index A numerical study followed 



in 



651 ] . where one peculiar two field model was investigated (backreaction was neglected) 



and indeed a slight enhancement was found. However, in the presence of many inflaton fields 



it was recently shown in 68j that resonances are generally suppressed: if inflatons dephase, 



which occurs quickly in most models even if the fields start to oscillate in unison, an increase 
of the matter field's effective mass results, because m^-faff) contains a term proportional to 
oc Yli V 9 ? • Since the matter field turns heavy and the oscillations of the matter field's effective 
mass are smeared out, it becomes more difficult to produce x particles. As a consequence, 
resonances are suppressed. This suppression is strong for large Af, as seen in tie case of Ag- 
nation where M ~ 100 [ssl; this numerical result was confirmed in In 68] preliminary 
results were reported, indicating that fewer fields would still allow for resonances, which 



65|]. 



appeared to be weaker than in the single field case, in contrast to the findings of 

Here, we provide a complete numerical simulation, including the expansion of the universe 
and backreaction effects, to settle this issue. It is evident from Fig. [1] A that resonances 
are indeed suppressed, even in the presence of only two fields. Thus, the initial hope that 
particle production would be enhanced in Cantor preheating is not true in general. The 
slight enhancement found in 65| might have been caused by the chosen coupling constants 
in conjunction with ignoring backreaction. Further, for as few as three fields, parametric 
resonance is heavily suppressed and disappears entirely for as few as five fields over the 
time-scale of our simulations (we cover over one hundred oscillations of the inflaton fields). 

The onset of parametric resonance shifts to later times if M increases from one to two, 
but stays about the same if M is further increased, Fig. [I] A. This effect could be caused by 
our chosen initial conditions, which do not coincide precisely with the breakdown of slow roll 
if more than one inflaton is involved. We cannot exclude particle production on longer time 
scales, but since the universe is expanding all along, causing the energy in the oscillating 
inflaton fields to redshift, it becomes more difficult, and at some point impossible, to start 
preheating. As a result, the universe remains dominated by the oscillating inflaton fields 
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until they decay perturbatively (i.e. mediated by gravity). 

The variance (IA22I) of the matter and inflaton fields is plotted in Fig. [2j A variance of 
order one indicates that a field is fragmented. We see that fields quickly become inhomo- 
geneous, and backreaction becomes important, once explosive particle production occurs. 
For M > 3 no fragmentation of the fields is evident, simply because particle production is 
absent in the first place. 

We conclude that the well known parametric resonance models based on four-leg inter- 
actions (gip 2 x 2 /2) are unlikely to work without fine tuning if more than two inflaton fields 
contribute to preheating. This conclusion is in line with the analytic arguments given in 



68|, |69(, but could be evaded by coupling each inflaton field to its own preheat matter field. 



However, such a construction seems somewhat artificial and ill motivated to us. 

In the absence of preheating the old theor y o f reheating (perturbative decay over 

many oscillations) still applies, as proposed in [68( . However, a delayed decay of the inflatons 
might pose problems if a field survives long enough to interfere with the thermal history 
of the universe. This problem was pointed out in [77|] within A/"-flation, and tachyonic 
preheating H[ was proposed as an alternative; this type of preheating should be less affected 
by dephasing, since the matter field's effective mass contains a term proportional to oc (fi, 
which is still negative half of the time. We investigate this interesting possibility in the next 
sections. Our main concern is that the incorporation of a four-leg interaction (g ^ 0) in 
addition to a three-leg interaction (a 7^ 0) could again spoil preheating. 

We would like to point out that a suppression of preheating might be desirable: the 
reheating temperature after preheating is generally high enough to violate bounds originating 
from preventing gravitinos from over-closing the universe or spoiling the success of primordial 
nucleosynthesis jlOo| . The maximum reheating temperature is usually taken to be T r h < 
W 9 GeV (g, 3 (see i.e. Q for a pedagogical review). This problem can be resolved via a 
second phase of reheating: if a long lived light scalar field dominates the energy density of 
the universe after the decay of the inflaton(s), i.e. because it redshifts like matter due to 
oscillations in a quadratic potential, it causes a second phase of reheating once it decays; 
the dilution of any particles present prior to this second phase of reheating alleviates the 
strict bounds on the reheating temperature 55 1. A concrete realization of this scenario is 

n n 

the decay of the saxion within the framework of F-theory [56( (see also [62J for related work 
withing G2-MSSM models). Given a multi-field inflationary scenario, the role of this late 
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decaying field may be played by one of the inflatons that did not decay successfully during 
preheating. 



B. Generic Preheating Case (a ^ 0, g ^ 0): no Suppression 

In addition to the four-leg interaction, let us turn on a three- leg interaction (a 7^ 0). 
Consequently, the square of the effective mass picks up a contribution proportional to 
a J2i L Pi + 9J2i i P'i- The first term oscillates around zero with a damped amplitude (due 
to dephasing) potentially leading to tachyonic preheating, while the second one resembles a 
positive offset that redshifts like matter. 

Due to the tachyonic term, we see an enhancement of preheating in Fig. [T]B as compared 
to panel A (a = 0) (compare also the variance in the left column of Fig. [3] with Fig. [2]); 
there is no discernible suppression of resonances if the number of fields is increased, but we 
still observe the slight delay in the onset of particle production if M is increased from one 
to two. 

We conclude that incorporating a tachyonic contribution leads to efficient preheating 
after multi-field inflation, even in the presence of sizable four-leg interactions. 



C. Tachyonic Preheating (er ^ 0, g small): Slight Enhancement 

Here, we investigate whether tachyonic preheating caused by a three-leg interaction 
(a 7^ 0) is affected by t he pr esence of several inflaton fields if the four-leg interactions 



are tuned down (g small) 101] . Since the matter field's effective mass is dominated by a 
term proportional to Ylii^Pi instead of Yli^fli we do n °t expect any suppression: the de- 
phasing of the fields leads to a damping of the overall amplitude of which remains 
negative half of the time. 

As evident from Fig. [I] C, an instability appears fast, leading to efficient preheating 
within a few oscillations of the inflaton fields. The variance is plotted in Fig. [3] (right 
column) . Interestingly, we observe a slight enhancement with an increase in M: for a single 
field, a few spurts of particle production, that is more than one oscillation of the inflaton, 
are needed to reach p x / p v of order one. But, for M > 3 this ratio approaches one following 
the first spurt of particle production and the fields fragment fast, as evident in Fig. [31 As a 
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consequence, the potential for the inflatons does not need to curve upwards again for <p < 0, 
but may simply flatten out. This enables preheating after inflation driven by fields with 
runaway potentials and could ease moduli trapping on the landscape [151 ]. 



III. DISCUSSION AND CONCLUSION 

Our lattice simulations show that preheating via four-leg interactions is strongly sup- 
pressed if more than two inflatons couple to the same matter field. The cause of this 
suppression is a dephasing of the inflaton fields that leads to an increase of the matter field's 
effective mass, while smearing out oscillations in it. One could imagine that each inflaton 
couples to a different preheat matter field, which in turn decay to standard model particles 
in the end. However, such a scenario seems fine tuned in the presence of several hundred 
inflatons, such as in jV-flation; even if only three inflatons, out of many, happen to be cou- 
pled to the same matter field, they would be unable to decay during preheating, and could 
dominate the energy of the universe; this might be useful, since their decay, i.e. mediated by 
gravity, could alleviate the gravitino problem. Fields are often of the same kind in current 
models of multi-field inflation, for example, they are all identified with axions in jV-f lation. 
Therefore, it seems rather unlikely for inflatons to have different decay channels. lOij ] 

We also considered tachyonic preheating where three-leg interactions (Yukawa couplings) 
lead to a tachyonic contribution to the matter field's effective mass, causing explosive parti- 
cle production even if the inflatons run out of phase. In the presence of additional four-leg 
interactions and more fields we find no suppression. Further, if four-leg interactions are 
suppressed, we observe a slight enhancement in the efficiency of tachyonic preheating with 
increasing M. Due to the employed tachyonic instability and the multiple fields involved, 
there might be additional observational signatures: the spectrum of gravitational waves pro- 
duced during preheating might differ and additional non-Gaussianities might be produced. 
We leave these interesting avenues for future studies. 

To conclude, tachyonic preheating is a viable mechanism to transfer energy to other 
scalar fields after multi-field inflation, opposite to standard parametric resonance scenarios; 
in the latter case, we expect some (or all) fields to survive long enough for the old theory of 
reheating to commence. 
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APPENDIX A: IMPLEMENTATION OF LATTICEEASY 

The current model, omitting the trivial Af = 1 case, consists of Af inflatons, tpi, each of 
which has 

2 Af + 2% - 3 2 2 , A -n 
m i = 2 (JV- l) m = ' ^ ' 

which describes a system that has an average square mass of 

Af 

AT 



nr " : 77 X>™ 2 ( A2 ) 



i=i 

N 



E<^+ 2 ;-^ (A3) 



1 2Af{N-l) 

m 1 (A4) 



A/' 2 + Af(A/'+i)-3A/'_ 2 



2Af(N - 1) 

= m 2 , (A5) 

as well as a lowest square mass of m\ = m 2 /2, and a greatest square mass of m 2 ^ = 3m 2 /2. 
The total potential energy is equally distributed among the Af fields and sums up to m 2 ip 2 ) /2, 
where tpo = 0.193 lo| . in line with single- field models. The resulting initial potential energy 

iAmV 2 (0) = -limV 2 , (A6) 

yield the initial conditions 



Vi{$) = \l^jW, (A7) 
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for the zero modes of the inflaton fields. We additionally assume that they are initially at 
rest, so that the homogeneous momentum vanishes, 



V?i(0) = 0. 



(A8) 



The inhomogeneous initial conditions are naturally defined in momentum space and na- 
tively calculated by LatticeEasy. The Fourier transform of any field in the model, fc, is 
parameterized by an amplitude and a phase 



/<(*) = /*(*) 



(A9) 



where we use the convention 



(A10) 



The amplitude of each mode is chosen from a Raliegh distribution 



za 



73 



93), 



V 



fi(k) = fi(k) 



-2a 2 u) k \fi(k)\ 



where 



and 



ujk = k 2 + a 2 m 



icS 



m 



ieff 



(All) 

(A12) 
(A13) 



The phase Qi is taken to be evenly distributed between and 2n. 
The derivative of the field is determined analytically, 

fi(k) = ±iu k fi - Hfi, 



(A14) 



where the ambiguity in sign comes from choosing either a right-moving or left-moving wave. 
Although each Fourier mode of (the real field) fa should match the Raliegh distribution, 
simply choosing modes from (1A1 1[) does not guarantee that the field and its derivative are 
real while simultaneously preserving isotropy. Namely, the modes of the field must obey 



f i (-k) = f*(k). 



(A15) 



Similarly, the modes of the derivatives must obey an analogous relation. Felder and Tkachev 
impose this condition by hand, while preserving the initial spectrum of the field, by choosing 
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two modes whose amplitudes are drawn from (IA11I) with a random phase, and fi^ik), 

and defining the initial conditions of the field to be 

(A16) 



Mk) = \l-[ki(k)+f h2 (k) 



while the initial derivatives of the fields are 



fi(k) = \j-iuJk \fi,i(k) - fi, 2 (k) ) - Hfi(k) . 



(A17) 



Lastly, we rescale time, space and fields so that the homogeneous fields and the fields' 
derivatives are of order one, resulting in our definition of a program time 



a program comoving distance 



and program fields 



r = mt . 



a 3 / 2 

fi,pr fi 



(A18) 



(A19) 



(A20) 



With the above rescaled fields, we define the average fields 

1 " 3 

and the variance 



(A21) 



3=1 



Var ( 



\ 



n 3 (^'P r ^P r ) 



(A22) 



3=1 



where f? is the field value of the i'th field in the j'th lattice cell (we use an n 3 = 128 3 



lattice so that j = 1 . . . n 3 ). 
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[95] This problem arises in the heavy moduli scenario 
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moduli-problem 
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13, 



57|, 



581 ] . a solution to the cosmological 



[96] One can also reheat via MSSM flat directions 



611 ]; see also 



6J. 



[97] In some concrete scenarios the mass distribution is known (i.e. in AA-flation [331]). and should 



be used during preheating (see i.e. 



[98] Note that three-leg interactions are generated if ipi = Vi does not coincide with the location 

where x becomes lightest according to four-leg interactions. 
[99] Without a compelling theoretical reason, vastly different coupling constants constitute fine 

tuning. 

[100] Often the gravitino is the lightest SUSY particle with a long lifetime, providing a candidate 
for dark matter. However, its decay and the resulting shower of energetic photons can cause 
the disintegration of light elements whose abundances, as predicted from nucleosynthesis, are 
in excellent agreement with observations (except Lithium, which is off by a factor of 2 to 5). 

[101] Setting g = leads to instabilities in our code. 
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82, 



since reheating is 



[102] The situation is different in models of staggered inflation 

dominated by a few long lived fields which can be quite different from the already decayed 
ones (i.e. they could lie within the MSSM). 



20 




FIG. 1: We plot the ratio p x / ' p^ of the total matter energy density (kinetic, potential and gradient) 
to the total energy density in the inflatons = YliPvi ( n °te X 2 ^liig^Pi + cnpi)/2 is not included 
in either pu, or p x ) for cases A, B and C (n 3 = 128 3 -lattice, L = 5/m, for coupling constants see 
Table H]). Preheating is successful if this ratio becomes of order one; a turbulent phase and an era of 
thermalization follow (not covered here). Before p x j p v ~ 1, the fields fragment and backreaction 
becomes important, see Fig. [2] and [3] for the variance of the fields. Panel A: Parametric resonance, 
and thus preheating, is strongly suppressed for J\f > 3. Panel B: a comparison with panel A 
reveals that resonances are enhanced due to tachyonic preheating, which is caused by ax 2 
an increase in the number of fields has no noticeable effect, except for a slight delay of particle 
production. Panel C: tachyonic preheating is slightly enhanced if the number of fields is increased. 
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FIG. 2: The variance fA22l) of the fields for case A (a = 0, g + 0, 



n" 



128 3 -lattice, L = 5/m), 



increasing the number of fields from ]\[ = 1 to 5. Fields fragment and backreaction becomes 
important once the variance is of order one. For M > 3 the fields remain essentially homogeneous, 
due to the suppression of resonances caused by dephasing fields; preheating is inefficient in these 
cases. 
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FIG. 3: The variance fA22l) of the fields for case B (a ^ 0, g / 0, n 3 = 128 3 -lattice, L = 5/m) and 
case C (g suppressed compared to case C), for M = 1,3 and 5. Fields fragment and backreaction 
becomes important once the variance is of order one. a ^ leads to tachyonic preheating, which 
is efficient even in the presence of a large g. 
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